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Introductory remarks 


Paper describes a numerical FE based aeroelastic-acoustics analysis 
methodology 

This is followed by solution details of a 2-D airfoil and results correlation with 
known theoretical solution 

Also presented the 3-D wing case and related results pertaining to vibration, 
steady and unsteady flow (CFD), aeroelastic and aeroelastic-acoustic 
simulations 

Further associated solution results are presented for a numerically simulated 
unsteady pressure data 

Also SPL results from a SOFIA flight sensor data is presented in some detail 

Provides a discussion on implementation of these techniques in an existing FE 
software suitable for solution of complex, practical problems 

Discussions and concluding remarks 
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Numerical Formulation 


Structural (FE) free vibration matrix equation solving for co and 9 



Mu + Ku = 0 


The aerodynamic data are next computed by solving the Navier-Stokes (FE) 
equation 
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• Vehicle equation of motion is then cast into the frequency domain 


Mq + Cq + kq + f a (t)+f,(t) = Q (3) I 


where 

m(=& t M@) is the generalized mass matrix and similarly K,C ; / a M is the 
aerodynamic load vector; fXt) being the generalized impulse force vector 









Numerical Formulation (cont’d) 


Generalized impulse force vector f,(t) 


Earlier, the CFD code analysis results were verified with flight test data 
- Hyper-X vehicle and flight data comparison 





Numerical Formulation (cont’d) 

• Equation (3) may then be cast in a state-space matrix form as 


/ oir?i o im o i r o 

0 I q - -M 'K -M 'C q ~ _-M A f a (t)~ r M A f,(t) 


or 


x s (t) = A sl x s {t)+B a f{t) (5) 


where 


A„ = 



■ = x s = 



and 


y s (t)=c s ,x s (t)+Df(t) 


in which C=I and D =0 
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Numerical Formulation (cont’d) 


• In the presence of sensors, for the most general aeroservoelastic case these 
equations are converted into zero order hold (ZOH) discrete time equivalent at 
the k-th step: 


x s (k + l)=G s x s {k)+H s f(k) 

( 8 ) 

yXk + \)=e s x s {k)+Dj{k) 

( 9 ) 


in which 



where A s and B s are A st and B st having been modified to include sensors 
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Numerical Formulation (cont’d) 


• Coupled aeroelastic (AE) model 




Numerical Formulation (cont’d) 


• Acoustic frequencies are obtained by performing FFT on computed unsteady 
aerodynamic pressures 

• Also the sound pressure level (SPL) for a specified node is computed by first 
fixing a time band, t and then performing the following calculation using n 
number of sampling points 


a) compute average pressure 


f av 8 =(Z;^)A < 11 >[ 


b) compute the root mean square of pressure 




c) compute the SPL 



SPL(c/B)= 201og 


10 


rms 


ref 


(13) 


where p ref = 20xl0~ 6 Pa for air 

r~9 







Input structural/aero data 


Steady-state 

Euler/Navier-Stokes solution 
(local time stepping) 


Modal analysis 
(p's, CD's 

Form G S H S C S D S 
(structural state-space matrices) 


For each of a number of a specified 

p 


fa=*JbA 


Impulse load f x input 
for t = t s thru t e 

Form G S H S C s D s 
(structural state-space matrices) 
Solve for q,q (k+1 step) 

Control law to enforce desired 
control surface motion 
modified q,q 

u = 0 q 
u-0 q 

If u is large, regrid aero 


Change velocity B.C. at solid 
boundary feed u, U 


Unsteady aerodynamic analysis 
Euler/Navier-Stokes solution 
(one-step global time stepping) 
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Numerical Example I, 2-D 


2-D subsonic NACA 0012 Airfoil 

- Flight Condition : Mach Number = 0.3, Angle of Attack = 5-degree 

In an effort to verify the solution accuracy: 

- Steady state : Correlation with (i) Smith-Hess panel method (ii) CFD solver 

- Unsteady analysis solution is compared with Wagner’s suddenly accelerated airfoil 
problem 
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Numerical Example I, 2-D (cont’d) 


• Steady state solutions 


1 . Coefficient of Pressure (Cp) 
distribution over solution domain 


2. Coefficient of Pressure (Cp) 
distribution comparison 
on the airfoil 




Numerical Example I, 2-D (cont’d) 


• Unsteady analysis solution 

- Unsteady Coefficient of Lift (Cl) history comparison 





Unsteady analysis solutions 
- Coefficient of Pressure (Cp) 





Numerical Example I, 2-D (cont’d) 


• Unsteady analysis solutions 
- Velocity 
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Numerical Example II, 3-D 


• 3-D cantilever wing with NACA 001 2 airfoil 

- Flight condition : Mach 0.3 and 0.6, Angle of Attack = 0 degree 


1 . Cantilever wing with 
aeroelastic solution domain 


2. Structural (a) and 

aerodynamic (b) surface grid 
of wing 





Numerical Example II, 3-D (cont’d) 


• Structural modes of cantilever wing 



(a) mode 1, (b) mode 2, 

3.521Hz., IB 14.819 Hz., IT 


(e) mode 5, (f) mode 6, 

41 .408 Hz., 3B 48.636 Hz., 2T 



(c) mode 3, (d) mode 4, 

1 8.81 9 Hz., 2B 38.026 Hz., in plane 


(g) mode 7, (h) mode 8, 

59.797 Hz., 1 B control 64.958 Hz., 4B 






• Aerodynamic steady-state pressure and Mach distribution (M = 2.0) 




generalized displacement 


Numerical Example II, 3-D (cont’d) 




Numerical Example II, 3-D (cont d) 


Layout of the acoustic computation sampling points 
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Numerical Example II, 3-D (cont d) 


Mach 0.6 


1. Aeroelastic unsteady 
pressure response 
(node, C25) 


2. Acoustic wave 
frequencies 
(node, C25) 


Unsteady pressure history for Mach 0.6 
( C-25, Control Surface ) 



0.8 1.0 1.2 
time (sec) 


Acoustic wave frequencies for Mach 0.6 
( C-25, Control Surface ) 
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Numerical Example II, 3-D (cont’d) 

• Sound Pressure Level (SPL) for a line along the cord length (20% from wing tip) 




Introductory remarks 


Paper describes a numerical FE based aeroelastic-acoustics analysis 
methodology 

This is followed by solution details of a 2-D airfoil and results correlation with 
known theoretical solution 

Also presented the 3-D wing case and related results pertaining to vibration, 
steady and unsteady flow (CFD), aeroelastic and aeroelastic-acoustic 
simulations 

Further associated solution results are presented for a numerically simulated 
unsteady pressure data 

Also SPL results from a SOFIA flight sensor data is presented in some detail 

Provides a discussion on implementation of these techniques in an existing FE 
software suitable for solution of complex, practical problems 

Discussions and concluding remarks 




Numerical Example III, simulated data 


100-Sample Window, STARS Code 

f P_VAL P_AVG P_RHS P_SPL | 



0 5000 100QO 15000 20000 25000 30000 15000 40000 45000 50000 55000 60000 


Samples 

1000-Sample Window, STARS Code 

— P_VAL — P_AVG P_RMS F_SPL | 



Samples 



delta SPL(dB) SPL(dB) 


Numerical Example III, simulated data (cont’d) 
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Post Analysis - STARS Code 



10000 20000 30000 40000 50000 G0000 

samples 




Amplitude 


Numerical Example III, simulated data (cont’d) 
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pressure 


Numerical Example IV, 

SOFIA measured unsteady pressure data 



(a) Measured sensor data 



SOFIA sensor measured unsteady pressure data 





SPL 


Numerical Example IV, 

SOFIA measured unsteady pressure data (cont’d) 




Numerical Example II, 3-D (additional) 


Mach 2.0 
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STARS capability 
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Concluding Remarks 


An integrated finite element based aeroelastic-acoustics analysis algorithm 
presented 

Its implementation in a code, suitable for large scale computation, is also 
presented in some detail 

Numerical verification example problem is demonstrated 

A 3-D wing problem is analyzed in detail, that demonstrate the FE codes 
capability to solve practical problems routinely 

Also presented a simulated example problem followed by a SOFIA flight 
measured data solution 







